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We present results of simulations of Two Color QCD using two flavors of Wilson quark in the 
fundamental representation, at non-zero quark chemical potential jj,, on an 8 3 x 16 lattice. Results 
for the quark number density, quark and gluon energy densities, and superfluid condensate are 
qualitatively distinct from the behaviour expected on the assumption that the dominant degrees of 
freedom are tightly bound scalar diquarks which Bose condense; rather the scaling with ji is more 
suggestive of a Fermi surface disrupted by a Cooper pair condensate. We also present evidence 
both for screening of the static potential, and color deconfinement, arising solely as a result of a 
non-zero quark density. 
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1. Preamble 

QCD with gauge group SU(2) and quark chemical potential fi ^ has enjoyed much recent 
interest [j|]. It lies in a class of models with long-ranged interactions (together with QCD with 
non-zero isospin chemical potential (ie. \L U = and vectorlike theories with adjoint quarks) in 
which the quark determinant detM(/i) is real, and hence amenable to study using standard LGT 
algorithms. 

In Two Color QCD, q and q fields live in equivalent representations of the color group SU(2), 
implying that chiral multiplets contain both qq mesons and qq baryons. In the chiral limit where 
m n -C m p , where p generically denotes any non-Goldstone hadron, it is possible to study the }i- 
dependence of the model systematically using chiral perturbation theory (#PT) [Q]. The key result 
is that the ground state has non-zero quark density n q > for all /I greater than some onset value 
Ho = \m n . At the same point a superfluid order parameter (qq) ^ develops, signalling spon- 
taneous breakdown of the original global U(l) baryon number symmetry. Since the transition is 
second order the matter, consisting of tightly-bound qq scalars, may be arbitrarily dilute, and in 
the limit jj, — > jx 0+ is a textbook Bose-Einstein Condensate (BEC). In this regime the quantitative 
predictions of #PT read: 

n q « fl(]X ~ Mo); (qq) « J I- (^j Hm (qq(fl)) = const. (1.1) 

This behaviour has indeed been confirmed by simulations with staggered fermions |g[ |], |||] 

This picture should be contrasted with another paradigm for superfluidity, namely BCS con- 
densation of weakly interacting quark Cooper pairs from opposite points of a Fermi surface with 
radius kp i=a Ep = pL. For such a system, we obtain n q simply by counting states within the Fermi 
sphere, and the condensate by assuming that only states within a layer of thickness A <C /I around 
the Fermi surface participate in the pairing: 

n q c<n 3 - ( qq )ocAn 2 . (1.2) 

The same rule of thumb predicts the quark energy density e q oc pc. 

We might expect the crossover from BEC to BCS for (ju — jU G ) ~> \ (m p — m n ), at which point 
the Goldstones are no longer distinguished hadrons, and the fermionic nature of the constituents 
comes into play. The chiral limit is thus unimportant, and for this reason we have chosen to re- 
visit Two Color QCD using Wilson fermions Pioneering simulations by the Hiroshima group 
appeared in [0|. As shown below, an advantage of the Wilson formulation is that simulations 
with Nf = 2 quark flavors are possible, ensuring the theory is asymptotically free and confining 
for ;ii = T = for all couplings, with a controllable continuum limit as j3 — > °°. This makes the 
Wilson LGT particularly suitable for a study of gluodynamics in the presence of a background 
baryon charge density, though it should be stressed that in contrast to 3-color QCD, there is here 
no physical distinction between fundamental and anti-fundamental charges. 

2. Formulation and Simulation 

The Nf = 2 quark action we want to simulate is as follows: 

S = xj/iMYi + ¥iMY2 -JViiCYs^ivZ + JVI{CY5)t:2V\, (2.1) 
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where flavor indices are written explicitly, and we have included scalar isoscalar diquark source 
terms, which will help to reduce IR fluctuations in the superfluid phase [gj]. The action coincides 
with that of ^ with the identification J = 2], J = 2j. The charge conjugation matrix C defined 
by Cy^C -1 = 7^ with 7^ = in Euclidean metric has the properties = C 1 = —C, [C, 75] = 0. 
The crucial Two Color identity is %2V ii {x)X2 = U*(x). The source strengths J and / are a priori 
independent and their signs are arbitrary for now. The matrix M is the textbook Wilson fermion 
operator 

Mxy(n) = 8 XJ - k£[(1 - 7 v )^ [/ v (x)S M+ v + (1 + 7v)e-^ Sv0 Ul(y)8 y ^y] (2.2) 

V 

with properties 

7 5 M(At)75 = M\-/i) (2.3) 

Cz 2 M{ix)C- l Z2 = M tr {-n) (2.4) 

(C7 5 )t 2 M( j u)(C7 5 )- 1 t 2 = M*(n) (2.5) 

Property ( |2~3| ) implies detM(/x) is real, but because M contains both hermitian and antihermitian 
non-constant components there is no proof that it is positive. This confirms the impossibility of 
simulating a single Wilson flavor without a Sign Problem, and also suggests that even for Nf even 
there is an ergodicity problem along the lines of that found for adjoint staggered fermions [H], where 
changing the sign of detM requires an eigenvalue to flow through the origin, which the Two-Step 
Multi-Bosonic algorithm can manage but the Hybrid Monte Carlo (HMC) algorithm not. 

Now, with the change of variables = — xj/tfCfy, = C^ X X2^{ and the dropping of the index 



from flavor 1, eq. (2.1 ) can be recast as: 



»,*)C£ ) ..'*)(?)=***. (2-6) 



-7"7 5 M(-n)J V0 

Note that the action is bilinear in the variables so that the Grassmann integral yields a factor 
det^ rather than a Pfaffian. Using the identity 

det( 7 J =detXdet(Z-WX- l Y) (2.7) 



W Z 

and properties ( |13| , [14"| ) we deduce 



deUf = det(M t ( i u)M( J u )+JJ). (2.8) 

Hence positivity of det^# requires the product J J to be real and positive, which translates into 
the requirement that the diquark source term be antihermitian [§J. Since no eigenvalue of ^# can 
vanish, the ergodicity problem is also cured. 



Now use (2.6) to write 



{M*(ji)M(jl) + \J\ 2 



JtJt= t»-WJ»WTV\ . (2 . 9) 

V M t (- J u)M(-At) + |7| 7 1 
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The off-diagonal terms can be shown to vanish if J = J* using (|2J|); moreover the same identity 
applied to the lower block yields 

deUf 1 ^ = [det(A^(ji)M(ji) +JJ)] 2 . (2.10) 

It is therefore possible to take the square root analytically, by using pseudofermion fields with 
weight (MtM+ |/| 2 ) -1 . This has the advantage of (a) requiring matrix/vector multiplications of 
the half the dimensionality, and (b) permitting a hamiltonian evaluation and hence the use of an 
exact HMC algorithm. It is equivalent to the even/odd partitioning step used for staggered fermion 
gauge theories, but is more transparent since all lattice sites are physically equivalent, making the 
force term easier to implement. The trick was used in [9], though because the staggered version 
still requires a Pfaffian rather than a determinant, an HMD algorithm was used in that case. 



3. Observables 

The most fundamental is the quark density, obtained as a derivative of the free energy: 

-K(Uft ~ l)e^U (x)<t> x+() + UYo + l)e-^(x-6)<!> x _ 6 ). (3.1) 

The apparently unphysical term — 2fc(y/(do — ju)i//) is irrelevant in the long wavelength limit. We 
have checked that at saturation hm^^Uq = 2NfN c . The quark energy density e q is approximated 
by the same expression with a relative minus sign between forward and backward timelike links. 

For the diquark condensate, it is convenient first to introduce rescaled source strengths {j,j} = 
fc -1 {J, J}, and then the orthogonal combinations j± = j ±j. We then have 

1 d In Jf K , _ — , ,„ „ x 

(??±) = -y-g— = 2^75<^T075^)- (3-2) 

Since the diquark condensate is not a component of a conserved current, its normalisation is to 
some extent arbitrary: we prefer to normalise it with the same factor of K as the quark density. 



4. Results 

We have studied an 8 3 x 16 lattice with parameters j8 = 1.7, K = 0.1780, using a standard 
Wilson plaquette action. Unfortunately our results are not directly comparable with those of [[7|] 
due to their use of an improved gauge action. Studies of the pL = string tension yield a lattice 
spacing a = 0.220fm, and the spectrum m n a = 0.800(2), m K /m p = 0.920(3) In simulations 
with n 7^ we have used a diquark source ja = 0.04. So far we have accumulated roughly 300 
HMC trajectories of mean length 0.5 for \xa G [0.3,0.9]. 

In Fig. 1 we plot quark observables (with a suitable subtraction so that e q (\L)= 0), and in Fig. 2 
the gluon energy density £ g = ^(tr(rj r — Qj) and Polyakov line, by far the noisiest observable. 
The behaviour of the superfluid condensate (qq) vs. jj, is qualitatively very different from the 
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Figure 1: Quark observables 



Figure 2: Gluon observables 
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Figure 3: Quark number density vs. [Xr 
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Figure 4: Diquark condensate vs. fi^ 



negative curvature seen in the j • —> limit in studies with staggered fermions[^, ^, g]. Moreover 



the magnitude of is much smaller, suggesting following (LI) that the effective f n is much 



smaller, and hence the hadron degrees of freedom more strongly interacting. 



This motivates us to replot the data following ( |1.2[ ), as shown in Figs. 3 & 4. There is some 
evidence to support the formation of a Fermi sphere; n q increases if anything more rapidly than 
/I 3 , although the impact of lattice artifacts has still to be determined. Similarly, the condensate is 
consistent with BCS behaviour being recovered in the j — ► limit, though data at smaller j will 
be required to confirm this. The most spectacular evidence for a Fermi sphere comes from the 
energy densities, plotted in Fig. 5. A significant fraction of the energy density, O(30%), comes 
from gluons, suggesting that the medium which forms is strongly-interacting. Finally in Fig. 6 we 
show the static quark potential V(r), showing clear evidence for screening for \ia > 0.4. 

5. Outlook 

The parameter set chosen for our initial study appear to place the Two Color quark medium in a 
qualitatively different regime from that described by #PT and studied numerically using staggered 
fermions. If our guess that the system has a Fermi surface disrupted by a diquark condensate is 
borne out by further analysis, then the system may have a lot more in common with high density 
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Figure 5: Quark and gluon energy densities vs. fi 4 
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Figure 6: Static potential V(r) for various ji 



QCD than originally thought. One exciting possibility is that the onset transition predicted by 
#PT at \i () a ~ 0.4 does not coincide with the deconfining transition, which the Polyakov loop data 
of Fig. 2 suggest may not occur until \ia as 0.6(1). Future study will also focus on the hadron 
spectrum, inspired by the evidence for decrease in mass for vector states reported in and a 
study of the gluon propagator. 
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